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Measuring spin of a supermassive black hole at the 
Galactic Centre — Implications for a unique spin 
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ABSTRACT 

Wc determine the spin of a supermassive black hole in the context of discscismology by 
comparing newly detected quasi-periodic oscillations (QPOs) of radio emission in the 
Galactic centre, Sagittarius A* (Sgr A*), as well as infrared and X-ray emissions with 
those of the Galactic black holes. We find that the spin parameters of black holes in Sgr 
A* and in Galactic X-ray sources have a unique value of » 0.44 which is smaller than 
the generally accepted value for supermassive black holes, suggesting evidence for the 
angular momentum extraction of black holes during the growth of supermassive black 
holes. Our results demonstrate that the spin parameter approaches the equilibrium 
value where spin-up via accretion is balanced by spin-down via the Blandford-Znajek 
mechanism regardless of its initial spin. We anticipate that measuring the spin of black 
holes by using QPOs will open a new window for exploring the evolution of black holes 
in the Universe. 
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1 INTRODUCTION 

The Galactic centre, Sagittarius A* (Sgr A*), is a com- 
pact source of radio, infrared, and X-ray emissions having 
variability in the range of a few tens of minutes to hours 
(Baganoff et al. 2001; Genzel et al. 2003; Yusef-Zadeh et 
al. 2006). These emissions seem to originate from a hot and 
low-density accreting gas plunging into a supermassive black 
hole (Yuan et al. 2004; Kato et al. 2009). A precise measure- 
ment of its mass and spin is a long-standing issue for astro- 
physics to investigate the mechanism of energy extraction 
from spinning black holes for astrophysical jet production 
as well as the evolution of supermassive black holes along 
the cosmic hitory (Bardeen 1970; Blandford & Znajek 1977; 
Wilson & Cobert 1995). Although the mass of Sgr A* has 
been constrained by using the stellar orbit method, a pre- 
cise measurement of its spin for the best-estimated mass has 
been poorly conducted. 

Recently, Miyoshi and colleagues have detected multi- 
ple quasi-periodic oscillations (QPOs) of radio emissions in 
Sgr A* (Miyoshi et al. in prep.), whose periods are close to 
the Keplerian period at the innermost stable circular orbit 
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of a supermassive black hole with mass 4 x 1O 6 M0. Because 
of the excellent spatial resolution of the Very Long Baseline 
Array (VLB A), the quasi-periodic radio emission certainly 
originates from within the central sub-mas scale, approxi- 
mately 100 r g around the central black hole at a distance 
of 7.6 kpc, where r g = GM/c 2 = 0.01 (M/1O 6 M ) AU is 
the gravitational radius (G, M, and c are the gravitational 
constant, the mass of black hole, and the speed of light, re- 
spectively). This is the first time that such multiple QPOs 
have been identified in the vicinity of a supermassive black 
hole. The spatial pattern of emission regions cannot be ex- 
plained by the Keplerian rotation of a single emitting body 
at a given radius. 

Four simultaneous QPOs (16.8, 22.2, 31.4, and 56.4 
min) are detected and the first three periods are identical 
to QPOs in the near infrared and X-ray observations during 
different observation epochs (see Table 1). Three identical 
periods in the different wavelength are stable at least for 
several years and the frequency ratio of last two periods is 
close to 3:2. Such a stable double peak QPO is a well-known 
feature for high-frequency QPOs (HF-QPOs) in Galactic X- 
ray sources (Remillard & McClintock 2006). The multiple 
periodicity and their coincidence between the different wave- 
lengths, and also the different observation epochs, indicates 
that the origin of QPOs in the Galactic centre is closely re- 
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lated to the dynamics of an accretion disc feeding the black 
hole. Therefore we measure the spin parameter of a black 
hole in Sgr A* by using the period of QPOs based on disc- 
seismology (e.g., Nowak & Wagoner 1993). 



2 METHOD AND MODEL 

One promising mechanism of generating multiple QPOs is 
a global disc oscillation excited by the resonance between 
geodesic modes of the disc (the so-called resonant disc oscil- 
lation model: Abramowicz & Kluzniak 2001; Kato & Fukuc 
2006; Kato et al. 2008). The resonant frequency is the com- 
bination among geodesic frequencies at the radius where the 
resonance occurs. When the resonance condition is specified, 
both the resonant frequency and the resonant radius is de- 
termined uniquely in terms of the black hole mass M and the 
spin parameter a* = Jc/GM 2 where J is the angular mo- 
mentum of the black hole. Therefore the metric of the black 
hole can be constrained by the frequency of the QPOs. 

Resonance may occur at a radius where the frequency 
ratio of the geodesic modes is a ratio of small integers and 
resonant response can either spontaneously grow or damp 
the oscillation itself (Abramowicz & Kluzniak 2001). One 
of the most prominent resonances is a mode-coupling be- 
tween acoustic waves and non-axisymmetric modes such as a 
warp in the disc, the so-called wave- warp resonance (Kato & 
Fukue 2006; Kato et al. 2008). For example, this resonance 
is excited at a radius r res where Qk = 2k. Here Qk and 
k are the Kepler frequency and the epicyclic frequency, re- 
spectively (see Fig. 1 of Kato & Fukue 2006 for the rela- 
tion between a* and r rea ). Qk and k at the resonant radius 
r rcs = r ICS /r s measured at infinity are expressed as 
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as derived by Okazaki et al. (1987). The resulting frequencies 
of QPOs are rnQn ± K and uiQk where m is the azimuthal 
mode number, and some lower mode oscillations related to 
such resonances are reported by numerical studies (Kato 
2004). 



3 RESULTS 

3.1 Unified model of QPOs 

Figure 1 shows the period of the observed QPO overlayed 
with lower mode (m = 1,2) resonant periods related to the 
wave-warp resonance as a function of the black hole mass 
ranging from a stellar mass black hole to a supermassive 
black hole (skipping over the intermediate mass region). 
QPOs in the Galactic centre are selected with regard to 
the multiple detection among different wavelengths (Table 
1). We found that such QPOs in Sgr A* detected at identi- 
cal frequencies are consistent with a mass-period relation for 
the spin parameter a, ~ 0.4 (see Fig. lb). At the same time, 
HF-QPOs in the Galactic X-ray sources agree well with the 




Figure 1. Observed QPO periods as a function of black hole 
mass, (a) HF-QPO periods of different sources are shown as 
crosses with horizontal bars indicating the range of black hole 
mass (Abramowicz & Kluzniak 2001; Homan et al. 2003; Orosz 
et al. 2004; Remillard & McClintock 2006; Kato et al. 2008). 
(b) QPO periods of Sgr A* in different energy bands are shown 
(see Table 1). The black hole mass is assumed to be (3.7 ± 1.5) X 
10 6 Mq (Schodel et al. 2002). Resonant oscillations for m = 1 and 
2 are shown as solid (Ok), dashed (Qk + K )> dotted (f^K — «), 
and gray solid (2Qk) lines. Note that 2£^k — n = Qk + K and 
2f2K + k are omitted for simplicity. Thin and thick lines indicate 
the periods for the spin parameter a* = 0.3 and 0.4, respectively. 



resonant periods for the same spin parameter within the er- 
ror of the estimated mass (Fig. la). Therefore we identify 
the three identical periods (16.8, 22.2, and 31.4 min) with 
resonant modes 2Qk, Qk + k, and Qk, respectively. 



3.2 Unique spin parameter 

Now we can determine the spin parameter of black holes 
by using the periods of QPOs corresponding to Qk- For in- 
stance, 31.4 min is used for Sgr A* and periods of lower HF- 
QPOs are used for the Galactic X-ray sources. Note that the 
frequency of single peak HF-QPOs are treated as Qk ■ In or- 
der to constrain the resultant spin parameter, the estimated 
mass of a supermassive black hole in Sgr A* is taken from 
recent measurements (Schodel et al. 2002; Ghez et al. 2008; 
Gillessen et al. 2009). Figure 2 shows spin parameters of all 
samples evaluated by using the discseismic measurement. All 
spin parameters are relatively small ($J 0.7) in comparison 
with the equilibrium value of spinning black holes (f» 0.95) 
predicted by a numerical study (Gammie et al. 2004) . When 
all samples are fitted by using a linear relation as a func- 
tion of the black hole mass, the spin parameter becomes 
larger than 1 for black holes with M ^ 1O 7 M0. Instead of 
a linear relation, we obtain a best-fit unique spin parameter 
a* = 0.44±0.08 , which is depicted by a gray shaded region, 
for la uncertainty by linear least square fitting. 



3.3 Evolution of BH spin and mass 

Next, we should ask why black holes have a unique spin pa- 
rameter in spite of the fact that their age as well as mass 
accretion history may vary in general. Actually, our results 
contradict recent studies that predict extremely spinning 
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Table 1. QPOs detected in Sgr A*. 



Obs. epoch (UT) 



Obs. band 



Period (min) 



Ref. # 



2003/06/15 - 16 
2004/09 

2002/10, 2004.08 

2007/04/04 

2007/07/22 

2004/03/08 09:30 - 16:30 



K-band 
1.60, 1.87, 1.90 (£im) 
2-10 (keV) 
L-band 
L-band 
43 (GHz) 



16.8 ±2, 28.0 
33 ±2 
22.2 
22.6 
45.4 



Gcnzel et al. 2003 
Yusef-Zadeh et al. 2006 
Belanger et al. 2006 
Hamaus et al. 2009 
Hamaus et al. 2009 



16.8 ± 1.4, 22.2 ± 1.4, 31.2 ± 1.5, 56.4 ± 6 Miyoshi et al. in prep. 
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Figure 2. Spin parameters measured by discseismic method. 
Crosses indicate spins for the Galactic X-ray sources whereas di- 
amonds indicate those for the Galactic centre in terms of black 
hole masses measured by the stellar orbits method (Schodel et 
al. 2002; Ghez et al. 2008; Gillessen et al. 2009). A gray shaded 
region indicate the best-fit spin parameter a* = 0.44 ± 0.08 for 
Icr uncertainty. 



black holes (Shapiro 2005; Volonteri et al. 2005). In order 
to test the feasibility of such a small unique spin parameter, 
we have to study the spin-up process by mass accretion and 
the spin-down process by the energy extraction as a result 
of the Blandford-Znajek mechanism, simultaneously. 

Figure 3 represents the equilibrium value of spin and 
also the time evolution of black holes surrounded by a rel- 
ativistic standard accretion disc (Novikov & Thorne 1973; 
Page & Thorne 1974; see also Kato et al. 2008), assum- 
ing given disc parameters such as the viscosity parameter a 
(Shakura & Sunyaev 1973), the magnetized parameter /3, the 
ratio of the gas pressure to the magnetic pressure, and the 
mass accretion rate m = M/A/edd normalized by the Ed- 
dington mass accretion rates A/edd = AnGM / ck, bb where K eB 
is the electron scattering opacity (e.g., Kato et al. 2008). In 
general, these parameters are not independent because mag- 
netohydrodynamic (MHD) turbulence in the disc is thought 
to be the source of viscosity and their values can only be 
examined numerically. For instance, we employ a — 0.01 on 



the basis of three-dimensional MHD simulations showing the 
total stress corresponds to a ~ 0.02 — 0.06 (Hawley 2000; 
Machida et al. 2000) for m < 1 and a « 0.01 for m ~ 1 
(Hirose et al. 2006). Recent MHD simulations also exhibit 
the natural emergence of large-scale magnetic fields (the so- 
called magnetic tower) at the inner region of an accretion 
disc (Kato et al. 2004). The formation of a magnetic tower is 
key to the extraction of the energy and angular momentum 
of a spinning black hole by the Blandford-Znajek mecha- 
nism and it has been suggested that the necessary condition 
for the energy and angular momentum extraction at the in- 
nermost region of an accretion disk is ft f» 1 (McKinney & 
Gammie 2004). 

The equations we solved in this study are the followings: 



d\nM M 



V 



dt 



M Mc 2 
V 



-rr = Ml ia - — 
at SZf 



(3) 
(4) 



where Af, ei n , and h u are the mass accretion rate, the spe- 
cific energy and the specific angular momentum at the inner 
edge of the accretion disc, respectively. The electromagnetic 
power loss V from the black hole is assumed to be that of 
the Blandford-Znajek mechanism: 



-£7f (^h — ^f) 



(5) 



where th is the radius of the event horizon and Qf and S7h 
are the angular velocity of the magnetic fields permeating 
the horizon and the angular velocity of the black hole, re- 
spectively (see Moderski & Sikora 1996; Beskin et al. 2003). 
The strength of magnetic fields B± permeating the event 
horizon is assumed to be regulated by the pressure of ac- 
cretion disc pdisc so that B\ = 87rpdisc//3- Note that the 
electromagnetic power loss is not negligible when f3 is less 
than the order of the unity. 

The relativistic standard accretion disc model provides 
a complete set of equations for describing the pressure of ac- 
cretion disc at the given radius as a function of the viscosity 
parameter a, the black hole mass m = M/Mq, the spin pa- 
rameter a„, and the mass accretion rate m. For a given m, 
the radiation pressure dominated region appears within the 
radius: 

fb = r b /r g 

= 36a 2/21 m 2/21 m 1G/21 B~ 1(i/21 V 2/21 H~ 10/21 n l<i/21 <> 



L Q lb/21 (6) 



where B, T>, H, and Q are the general relativistic correction 
factors (Page & Thorne 1974). To summarize, the pressure 
of the accretion disc can be described as follows: 
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and 



Prad 
Pgas 



f < fb, 
f > fb, 



Prad 



= 1.4 x 10 16 (am)' 1 Tli dyne cm" 



-9/10 ■ 4/5t-, j 

m VC2 dyne cm 



(7) 

(8) 
(9) 



Pgas = 3.0 x 10 (am) 

where TZi = f^Zr 2 ©"^ and TZ 2 
p-si/aog-u/Sp-o/io^-i/agVB are the radial depen _ 

dence including the general relativistic correction factors 
at the Boyer-Lindquist coordinated radius f — c?r/GM. 
The radius for evaluating the strength of magnetic field is 
asssumed to be fo = 1.3f ms where f ma is the marginally 
stable circular orbit (Bardeen et al. 1972): 



fms = 3 + z 2 - {(3 - zi) (3 + zi + 2z 2 )} 1/2 , 
where 

^ = l+(l-a 2 ) 1/3 [(l + a ,) 1 /3 + (1 _ at) i/3] 



Z2 = (3a* + Z\ I 



(10) 

(11) 
(12) 



Finally, we rewrite the equation (0 & (© by using the 
normalized variables as 



din m 

dt 
da* 
~~dt ~ 



TEDD 
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(e irj 



t,bzJ 
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Am* 



(13) 
(14) 



where symbols are the Eddington time tedd = M/Medd, 
the specific energy input e"i n = e- ln /c 2 , the efficiency of the 
Blandford-Znajek mechanism tjbz = 'Pbz/A^eddC 2 , the spe- 
cific angular momentum input k n — ck u /GM, the horizon 

radius fu = c 2 rn/GM = 1+ f 1 — a 2 ) 1 2 , and k — Qp/ilu = 
1/2 for the maximum efficiency of the Blandford & Znajek 
mechanism. Here we assume that the inner boundary is at 
the marginally stable circular orbit and both the energy and 
the angular momentum of accreting matter at the boundary 
are advected into the black hole. The specific energy and the 
specific angular momentum at the boundary are: 



Gin — G m s — 



1 - 



3f n 



2V3 1 



2a, 



(15) 



(16) 



We numerically integrated equations (|13[) & (|14[) with given 
initial parameters and track the evolution of black hole mass 
and spin. We also determined the equilibrium spin for m = 
10, 10 6 , 10 s by solving da,/dt = in the equation JUJ by 
using bisection method. 

Figure 3a shows the equilibrium value of spin as a func- 
tion of am. The equilibrium spin becomes larger when either 
a or m becomes larger. The best-fit spin parameter deter- 
mined by the discseismic method corresponds to an equilib- 
rium value of m as 1. Figure 3b shows the time evolution of 
spin, where the spin parameter of each model converges to a 
unique value regardless of the initial one. When the mass ac- 
cretion rate is regulated by the Eddington value (m — 1), the 
spin converges to the equilibrium value m 0.55 for stellar- 
mass black holes within the order of 10 8 years and then 
slowly approaches the equilibrium value « 0.4 for massive 
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Figure 3. Time evolution of the black hole surrounded by the 
standard accretion disc with the suitable disc parameter a = 0.01 
and p = 1 for different mass accretion rates (rh = 0.1, 1.0, and 
10 denoted by a green, red, and blue line, respectively), (a) is the 
equilibrium spin parameter in terms of am for given black hole 
masses (M = 10, 10 6 ,10 8 M Q denoted by a black sold, dashed, 
and dotted line, respectively), (b) is the time evolution of spin 
parameter for the initial black hole mass Mo = IOMq with 
different initial spin parameters (a* = 0.0, 0.75, and 0.95 de- 
noted by a solid, dashed, and dotted line, respectively). A gray 
shaded region in (a) and (b) indicates the best-fit spin parame- 
ter a* = 0.44 ± 0.08 determined by the discseismic measurement, 
(c) is the time evolution of mass ratio M/Mq with different mass 
accretion rate and initial spin parameters. The curves are almost 
independent of the initial spin parameters. 



black holes. When m = 0.1, the spin converges to a value 
Rs 0.5 within the hubble time, but never actually reaches the 
equilibrium spin. Therefore the resultant spin is consistent 
with the small unique spin w 0.44 when the mass accre- 
tion rate is regulated by the Eddington value rh ~ 1 with 
the appropriate disc parameters. On the other hand, when 
the accretion disc is somehow in a super-critically accreting 
phase, with m — 10, the spin converges to the equilibrium 
value of ps 0.96 within the order of 10 7 years. Although 
the equilibrium spin of the super-critical accretion phase is 
larger than the unique value, it could approach to this value 
during the subsequent sub-critical accretion phase in less 
than 10 9 years. The evolution of the black hole mass is not 
affected by the initial spin parameter (see Fig. 3c). Note 
that the final mass becomes 10 6 times larger than the initial 
mass for m > 1. 



4 CONCLUSIONS 

It has been suggested that the supermassive black hole in 
the Galactic centre used to be in the nearly critical mass 
accretion phase for more than the order of 10 8 years. A pos- 
sible explanation for such a large mass accretion history is 
the massive star formations in the proximity of the Galactic 
centre region. During the critical accretion phase, the spin 
reaches the unique value and the mass becomes ~ 10 6 Mq, 
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which is then maintained during the subsequent low accre- 
tion rate phase. Note that stochastic mass accretion history 
may also help to create the moderately spinning massive 
black hole (King & Pringle 2006). Similarly, black holes in 
Galactic X-ray sources have been in the nearly critical accre- 
tion rate phase for order 10 s years as well, suggesting their 
companion stars should be low-mass stars. Because they 
have reached the quasi-equilibrium state, the limit-cycle ac- 
tivities and also the emergence of jets does not alter their 
spin evolution. Thus, we conclude that the spin parameter of 
a supermassive black hole in the Galactic centre has a unique 
value of a* — 0.44±0.08. Conversely, the mass of a black hole 
consistent with the unique spin is M= (4.2±O.4)xlO 6 M . 

Without detecting the event horizon, we have con- 
strained the mass and spin of the supermassive black hole 
at the Galactic centre. The method we used here depends 
entirely on geodesic frequencies that are independent of the 
distance and viewing angle of a black hole. Once the unique 
spin parameter of the black hole in the Galactic centre has 
been confirmed by detection of the event horizon in the fu- 
ture observations (e.g., Takahashi 2004), studies of QPOs 
in other galaxies will open a new window to survey the 
growth history of massive black holes (Markowitz et al. 2007; 
Gierliriski et al. 2008). 
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